Annihilation of singlet fermionic dark matter into two photons 
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1 Introduction 



Although there exist numerous cosmological and astrophysical indications which confirm 
the presence of the non-baryonic dark matter, the nature of this substance remains 
unknown [1,2]. One of the most important motivations for the extension of the standard 
model (SM) is to introduce an appropriate candidate for dark matter. Among the 
various candidates proposed by the theories beyond the SM, weakly interacting massive 
particles (WIMP) are the most popular one. For instance supersymmetry models with 
R parity [1, 2], the extra dimensional models with conserved Kaluza-Klein (KK) parity 
[3], the T-parity conserved little Higgs model [4], and so on provide a WIMP dark matter 
candidate. Meanwhile, these models introduce others degrees of freedom besides of the 
dark matter which have not been confirmed experimentally until now. Of course, in 
addition to the dark matter there exist others motivations such as neutrino oscillation, 
baryon asymmetry, hierarchy problem, and so on for considering such models. However, 
maybe it seems natural to consider the most economical model for the explanation of 
these anomalies. For instance, if we want to explain dark matter anomaly, adding a 
singlet scalar to the SM content to play the dark matter role is the most minimal model. 
Several authors have studied the singlet scalar dark matter [5]. Alternatively, adding a 
singlet fermion together with a singlet Higgs to the SM content leads to a renormalizable 
theory for explanation of the dark matter [6]. Losing the renormalizablity of the theory 
we can have more minimal model with singlet fermionic dark matter (SFDM) [7, 8]. In 
all above scenarios, the weak interactions of WIMPs are the main key for explaining the 
thermal production of them in the early universe (for review see [1, 2]). Additionally, 
these weak interactions can provide an opportunity for search dark matter candidates 
through their production in high energy accelerators [9], their direct detection [10], 
and their indirect detection i.e. astrophysical observations of the annihilation or decay 
products of the dark matter in our galaxy or beyond. 

In fact, the dark matter anomaly is a mutual problem between particle physics and 
astrophysics; on one hand, we need to know particle physical properties of dark matter 
to attribute a galaxy distribution to it and on the other hand we need to know the 
distribution of the dark matter to construe a direct or indirect detected signals as dark 
matter signals. 

In the WIMP scenario, the weak interaction of dark matter would produce observable 
SM particles, such as charged anti-matter particles, photons and neutrinos. Among 
these, neutrinos and photons have advantage in comparison to others because they 
keep their source information during the streaming. Moreover, the very small cross 
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sections of the neutrinos make their flux very difficult to detect. Therefore, the gamma 
ray signature of the dark matter has been investigated extensively (For review see [11] 
and references therein). Since the continuum gamma ray emission from dark matter 
annihilation could be confused with astrophysical backgrounds, e.g. emission from 
galactic cosmic rays or millisecond pulsars, the study of the monochromatic gamma ray 
is important. Monochromatic gamma ray signatures have been studied for some dark 
matter candidates in the literature [12]. In this paper we investigate monochromatic 
gamma ray in the context of the SFDM model. 

The most minimal model, which includes a singlet fermion to serve as dark matter, 
is discussed in [7]. In this model the leading interaction between dark matter and 
the SM particle is given by the dimension-five term Hipip. Here H, ip and A are 
the SM Higgs doublet, the dark matter fermion and a new physics scale, respectively. 
Moreover, it is possible to explain the dark matter production in the early universe 
by a renormalizable theory if the Higgs sector of the SM is extended by a new singlet 
Higgs [6]. Also in the noncommutative space-time we can have a singlet fermionic dark 
matter which is involved in the noncommutative induced U{1) gauge theory [8]. There 
are also some others phenomenological studies on this model [13]. 

In this paper, we reobtain the thermal average of the annihilation of SFDM into the 
SM particles. Then, demanding the correct relic abundance we derive the coupling of 
SFDM with the singlet Higgs versus dark matter mass. For this calculation we use a 
set of Higgs masses so that the recent CMS and ATLAS results are respected ^. Then 
we obtain the cross section of the annihilation of SFDM pairs into two photons. The 
thermal average of this cross section is calculated throughout a favorite WIMP mass 
interval including resonance region. We depict the SFDM annihilation cross section 
into two photons versus dark matter mass, m, using consistent with relic abundance 
condition. We have organized the paper as follows: In the next section, we review the 
renormalizable extension of the SM by a singlet fermion which can be served as a cold 
dark matter and derive the Qs consistent with the relic density abundance. In Section 3, 
we calculate the two photons annihilation of the singlet fermion pairs and its thermally 
averaged cross section. We summarize our result in the last section. 

■^http: / /cms. web. cern.ch/news/observation- new-particle- niass-125-gev 
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2 Renormalizable Extension of the Standard Model 
by a Singlet Fermion 

It is possible to consider a singlet Dirac fermion beyond the SM which can serve as cold 
dark matter. The lepton number and baryon number of this new fermion are taken to 
be zero. Hence, there is not any mixing between the new fermion and the SM fermions. 
If a singlet Higgs is added to this model, the interactions between the singlet fermion 
and SM particles can be explained by a renormalizable manner. Explicitly, the singlet 
fermion and singlet Higgs are coupled as follows: 

>Cint = gsi^i^S, (1) 

and the singlet Higgs can be coupled with the SM Higgs through the following terms: 

£sH = -XiH^HS - X2H^HS^. (2) 

Therefore, in addition to the above terms, we need to add the following terms for singlet 
fermion and singlet Higgs to the SM Lagrangian: 

= - m^)ip, (3) 

£s = l{d,S)id^S) - ^S^ - ^S' - (4) 

The scalar potential given in (2) and (4) together with the SM potential, -^i^WH + 
\o{WHY, lead to the vacuum expectation values 

and {S) = xq (5) 

for the SM Higgs doublet which gives rise to the electroweak symmetry breaking and 
for the singlet scalar Higgs, respectively. These quantities are written with respect to 
the parameters introduced in the Lagrangian by the following relations derived from 




the extremum conditions, |^ 



{H)=^ - and 1^ \ {s)=xo - 0: 



/i = Aofo + (^1 + >'2Xo)Xo, 

ml = -yxo - y^o - ^ - >^2vl (6) 

We define the scalar fields h and s as the deviation from the vacuum expectation values 
of H and S, respectively. As it is obvious these fields are not mass eigenstates. The 
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various components of the mass matrix are written as follows: 
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h=s=0 



h=s=0 



dhds 

The mass eigenstates hi and /12 are 



A3 ^ 
= (Ai4 



A 



4 2 



Ait| 
3 ""'^ 2x0 ^ 

2A2Xo)t'o- 



/ii = sin(6')s + cos{9)h, 
h2 = cos{6)s — sm{6)h, 



where the mixing angle 6 is defined by 

tan(^) E 



with y 



2m? 



. The Higgs boson masses mi and m2 are given by 



m?2 



(7) 



(9) 



(10) 



This extension of the SM introduces some new parameters. In the sense of experiment, 
there exist allowed regions in the parameters space where it is possible to explain the 
relic abundance of the singlet fermions so that they serve as cold dark matter [6]. 
Before the study of monochromatic gamma ray, we give a reanalysis of the relevant 
parameters. At the tree level of perturbation, the annihilation of the singlet fermions 
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into the standard model particles is given by [6] 
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where F^- is the decay width of hi for z = 1,2, and o"h,-/i^. , ahihjhk ^"^^ the annihilation 
cross sections of tptp into /ij/ij or hihjh^ with A; = 1,2. Now we can compute the 
thermal average of the cross section over s, {crv), which is given by [17] 



{av) 



Sm^TpKlim/Tp] 



r dsais){s-Am')^sKJ^), 

Jim? \^F/ 



(12) 



where Kn{x) is the modified Bessel function of order n. The freeze-out temperature Tp 
is estimated by the following iterative equation: 



Xp = In 



m 



/45M 



pi 



2%^ y 2g^XF 



'av) 



(13) 



where xp = m/Tp and Mpi = 1.22 x lO^'^^GeV is the Plank mass and = 91.5 [18]. 

We consider the SFDM with mass below 200 GeV which is the range of Fermi-Lat 
data in hand [14]. Moreover, we take the SM Higgs mass 125 GeV which is consistent 
with the new reported value by ATLAS and CMS, and 500 GeV for the other Higgs. 
Therefore, the annihilation of SFDMs into Higgs particles is not possible. There are a 
couple of free parameters in the thermally averaged annihilation cross section, in ad- 
dition to the dark matter mass and Higgs masses: the coupling between the singlet 
fermion and singlet Higgs, gs, the Higgses mixing parameter 6, and the coupling con- 
stants between singlet Higgs and the SM Higgs, Ai and A2. The dependency of Ai and 
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Figure 1: The figure shows the dependency of Qs to m with respect to the relic density 
condition. 

A2 occurs through the decay width of the Higgs messenger to the lighter Higgs in de- 
nominator, so in our desired range of energy these coupling constants are not relevant. 
This decay rate is also dependent on the mixing angle which is fixed in our calculation. 
Demanding the correct relic abundance for the dark matter we derive Qs in terms of m. 
Since in our calculations the cross sections directly depend on 5fssin(26'), we have fixed 
6 to conclude Qs (One should note that the effect of the appearance 6 in the Higgs decay 
width in the denominator is not very significant). Fig. 1 illustrates Qg vs dark matter 
mass with sin^ = 0.02, for instance. It is clear that, if one considers a stronger mixing 
then the values of gs should be smaller. It is notable that, the valley represents the first 
resonance region where the dark matter mass is about half of the standard model Higgs 
mass. In addition, the decreasing region begins at about 170 GeV where the threshold 
of top quark production is. 



3 The ipip — > 77 cross section 

Now, according to the model presented in the previous section, we derive the cross 
section for the annihilation of two dark matter particles into gamma ray pairs. Although 
our dark matter cannot directly couple to the SM particles, the annihilation can occur 
by a Higgs particle through s-channel. In addition, since photon is a massless particle, it 
does not couple to the Higgs boson directly. Nevertheless, H77 vertex, can be generated 
at the quantum level with loops involving massive particles, such as bosons and 
charged fermions. The dominant Feynman diagrams contributing to this process are 
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shown in Fig. 2. 




Figure 2: The dominant Feynman diagrams for the annihilation of two singlet fermionic 
dark matter particles into monochromatic gamma ray lines. 




Figure 3: Annihilation cross section into gamma ray lines as a function of dark matter mass 
for m/jj = 125 GeV, = 500 GeV, and gs coming from Fig. 1. Dashed line denotes the 
Fermi-Lat bound [14]. 

Since the final products are massless, we can write down the cross section as follows: 

spins 

where y/s is the energy in the center of mass frame, and 

j=12 '''"'hj'-j 

where 5*1^ and 5*2^ denote sin 6' and cos 9, respectively, and Mh.^^^ is the amplitude for 
the decay of a Higgs into two photons. Here, T j is the decay width of j'th Higgs. One 
can write Mh^^^^ as follows [15, 16]: 

yi-jS^eyags 2 
Mu,^,, [3(-) F, + Fwi 
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where 

Ft = -2r[l + (l-r)/(r)], 
Fw = 2 + 3r + 3T{2-T)f{r). 
Here, r = Arrii'^/s with i = t,W and 

sin^^A/l/rj , for r > 1 

-i(lni^^-^vr)' for r< 1. 

Therefore, we obtain the cross section as follows: 

spins V ^ ^ 

(14) 

Thermal average of this cross section is also obtained using (12) and (13). We have 
illustrated the phenomenological aspects of this calculation via the Fig. 3. The SFDM 
can be served as cold dark matter provided that the free parameters are chosen such that 
the relic abundance condition is respected. Therefore, this figure shows cross section 
versus the dark matter mass while we use Qs from Fig. 1. For comparison, we have 
also shown the upper limits obtained by Fermi-Lat [14] for Navarro- Frenk- White (NFW) 
density profile. It is noticeable that most of the region (except the resonance) is actually 
excluded by direct detection experiments, following the results of Figure 6 in Ref. [6] . 
However we should note that, for masses below about 70 GeV (and above about 150 
GeV), in particular the resonance regions where are not excluded by direct detection 
experiment, the cross section is comparable with the Fermi-Lat data. Therefore, these 
values may be explored by likely future more precise experiments. 




4 Summary and discussion 

In this paper we have considered an extension of the SM which includes a singlet 
fermion as a cold dark matter and a real singlet scalar as a messenger to explain the 
corresponding thermal relic density [6]. Demanding correct relic abundance, we have 
derived Qs and plotted it in Fig. 1 versus dark matter mass. This figure is consistent 
with figure 2 in ref. [6]. Only we have fixed the Higgs mixing angle to obtain a definite 
value for Qs- Then the cross section of the annihilation of the singlet fermion into two 
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monochromatic photons has been obtained. The thermal averaged of the annihilation 
of this cross section is calculated for a wide range of energies including resonance. We 
have depicted the thermally averaged of the SFDM annihilation cross section into two 
photons for a reasonable set of the Higgs masses in Fig. 3. This figure allows one to 
compare our result with the Fermi-Lat observations. The most of the region (except 
the resonance) is excluded by direct detection experiments (see figure 6 in ref. [6]). 
However, the Fermi-Lat data can be an alternative way for exploring the nature of dark 
matter. As one can see from Fig. (3), in some region (mass below 70 GeV and above 150 
GeV which include the resonance regions), the two-photon annihilation cross section of 
SFDM is comparable to the resent Fermi-Lat data and potentially can be explored by 
likely future more precise experiments. 
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